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Abstract 

We consider a Gaussian multiple-input multiple-output (MIMO) wiretap channel model, where there 
exists a transmitter, a legitimate receiver and an eavesdropper, each node equipped with multiple antennas. 
We study the problem of finding the optimal input covariance matrix that achieves secrecy capacity subject 
to a power constraint, which leads to a non-convex optimization problem that is in general difficult to 
solve. Existing results for this problem address the case in which the transmitter and the legitimate 
receiver have two antennas each and the eavesdropper has one antenna. For the general cases, it has been 
shown that the optimal input covariance matrix has low rank when the difference between the Grams of 
the eavesdropper and the legitimate receiver channel matrices is indefinite or semi-definite, while it may 
have low rank or full rank when the difference is positive definite. In this paper, the aforementioned non- 
convex optimization problem is investigated. In particular for the multiple-input single-output (MISO) 
wiretap channel, the optimal input covariance matrix is obtained in closed form. For general cases, we 
derive the necessary conditions for the optimal input covariance matrix consisting of a set of equations. 
For the case in which the transmitter has two antennas, the derived necessary conditions can result in 
a closed form solution; For the case in which the difference between the Grams is indefinite and has 
all negative eigenvalues except one positive eigenvalue, the optimal input covariance matrix has rank 
one and can be obtained in closed form; For other cases, the solution is proved to be a fixed point of 
a mapping from a convex set to itself and an iterative procedure is provided to search for it. Numerical 
results are presented to illustrate the proposed theoretical findings. 

Index Terms 
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Secrecy capacity, Gaussian MIMO wiretap channel, transmitter optimization, physical layer based 
security. 

I. Introduction 

Wireless physical (PHY) layer based security from a information-theoretic point of view has received 
considerable attention recently, e.g., Q-Hl, and the comprehensive overview in [6i|. Wireless PHY layer 
based security approaches exploit the physical characteristics of the wireless channel to enhance the 
security of communication systems. The wiretap channel, first introduced and studied by Wyner [Tj, 
is the most basic physical layer model that captures the problem of communication security. Wyner 
showed that when an eavesdropper's channel is a degraded version of the main channel, the source and 
destination can achieve a positive perfect information rate (secrecy rate). The maximal rate of secrecy 
rate from the source to the destination is defined as the secrecy capacity and for the degraded wiretap 
channel is given as the largest between zero and the difference between the capacity at the legitimate 
receiver and the capacity at the eavesdropper. The Gaussian wiretap channel, in which the outputs at 
the legitimate receiver and at the eavesdropper are corrupted by additive white Gaussian noise (AWGN), 
was studied in |8]. Along the the same line, the Gaussian MIMO wiretap channel was investigated and 
the secrecy capacity of the MIMO wiretap channel was established in terms of an optimization problem 
over all possible input covariance matrices ||9l, ifTOl . In |[9|, ifTOl . the Gaussian MIMO wiretap channel 
model was given as = HjX + rij, i = 1, 2 where iij is AWGN with zero mean and covariance cr^I, 
and the power constraint Tr(Ri.) < P was used, where R^, is the input covariance matrix. An alternative 
expression of secrecy capacity was derived in ifTTIl . |[T2|| for another Gaussian MIMO wiretap channel 
model, i.e., = x + Vj, i = 1, 2 where Vj is additive Gaussian noise (AGN) with zero mean and 
invertible covariance Wj. Further, in |[T2l . the power covariance constraint ■< S was used where S is 
a given matrix and R^; ■< S denotes that S — R^; is positive semi-definite, which allowed for the secrecy 
capacity to be obtained in closed form rather than as a solution to an optimization problem. 

For the former Gaussian MIMO wiretap channel model, i.e., yj = HjX + rij and the power constraint 
Tr(Rx.) < P, finding the optimal input covariance matrix that achieves secrecy capacity leads to a non- 
convex optimization problem. This problem is in general difficult to solve. The solution of a special case 
in which the transmitter and the legitimate receiver each has two antennas and the eavesdropper has one 
antenna was given in |13]. In |14], it was pointed out that the optimal input covariance matrix has low 
rank when the difference between the Grams of the eavesdropper and the legitimate receiver channel 
matrices is indefinite or semi-definite, based on the assumption that the Grams both have full rank. 
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In this paper, we investigate the aforementioned non-convex optimization problem. In particular, for 
the multiple-input single-output (MISO) wiretap channel, we obtain the optimal input covariance matrix 
in closed form. For general MIMO case, we derive the necessary conditions for the optimal solution 
consisting of a set of equations. Those conditions result in a closed form solution for ut = 2. For the 
more general case in which the difference of the Grams of the eavesdropper and the legitimate receiver 
channel matrices is indefinite and has all negative eigenvalues except one positive eigenvalue, we prove 
that the optimal input covariance matrix has rank one and can be obtained in closed form. Otherwise, we 
prove that the solution is a fixed point of a mapping from a convex set to itself and provide an iterative 
procedure to search for it. 

The remainder of this paper is organized as follows. The mathematical model is introduced in In 
^nijthe optimal input covariance matrix is obtained in closed form for Gaussian MISO wiretap channel. 
In ^ irvl we derive the necessary conditions for the optimal solution consisting of a set of equations for 
the Gaussian MIMO wiretap channel. In ^V) we obtain a closed form solution for the case in which the 
transmitter has two antennas. In WH for the case in which the difference of the Grams is indefinite and 
has all negative eigenvalues except one positive eigenvalue, we prove that the optimal input covariance 
matrix has rank one and can be obtained in closed form. Numerical results in Willi illustrate the proposed 
algorithm. Finally, ^IXl gives a brief conclusion. 

A. Notation 

Upper case and lower case bold symbols denote matrices and vectors, respectively. Superscripts *, 
T and f denote respectively conjugate, transposition and conjugate transposition. det(A) and Tr(A) 
denote the determinant and trace of matrix A, respectively. Amax(A) denotes the largest eigenvalue of 
A. A ^ and A ;^ mean that A is a Hermitian positive semi-definite and positive definite matrix, 
respectively. A ^ B denotes that A — B is a positive semi-definite matrix. A ^ denotes that matrix A 
is not positive definite. rank(A) denotes the rank of matrix A. diag(v) denotes a diagonal matrix with 
diagonal entries consisting of the elements of v. ||a|| denotes Euclidean norm of vector a, while ||A|| 
denotes Frobenius norm of matrix A. I„ denotes the identity matrix of order n (the subscript is dropped 
when the dimension is obvious), a ^ h means that a goes to b. Given a matrix A, the matrix A^A is 
called a Gram. 
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II. System Model and Formulations 

Consider a MIMO wiretap channel where the transmitter is equipped with ut antennas, while the 
legitimate receiver and an eavesdropper have n/j and ue antennas, respectively. The received signals at 
the legitimate receiver and the eavesdropper are respectively given by 

YR = Hi?x + UR, YE = H£;x + n^; (1) 

where {nR x ny), H^; {ue x ut) are respectively channel matrices between the transmitter and 
legitimate receiver, and between the transmitter and eavesdropper; x is the x 1 transmitted signal 
vector with zero mean and nx x ut covariance matrix ^ 0; and n^; are circular Gaussian 
noise vectors with zero mean and covariance matrices cr'^ina ^^'^ <^^InB' respectively. We assume the 
power constraint is P, namely, Tr(Rj:) = P. It is easy to verify that the problem under Tr(R2,) < P 
is equivalent to that under Tr(R^) = P. We consider the scenario in which the transmitter has perfect 
short-term channel state information (CSI). 
The secrecy capacity is defined as ifTOl 

Cs = max Cs(Rx) (2) 

R,hO,Tr{R,)=P 

where 

C,(R,) = log det(I„, + URK^U^a^) - log det(I„, + Hi^R^H^a^) (3) 
is the secrecy rate. 

The transmitter optimization problem is to determine Ra; that maximizes the secrecy rate, i.e., achieves 
secrecy capacity. The optimization makes sense when the secrecy capacity is positive. We assume 
hJjH/j — H^H^; 7^ since otherwise Cs(Rx) = 0. Whether Cg > depends on the difference between 
the Grams of the legitimate receiver and eavesdropper channel matrices, i.e., hJjH^ — H^H^;. The 
following lemma provides the conditions to maintain Cs > 0. 

Lemma 1: The sufficient and necessary condition for Cg > is: hJjH/j — H^H^; is positive semi- 
definite or indefinite. 

Please see Appendix |A] for details. 

Let us assume that the channel matrices H/j and H^; have been normalized so that Tr(HjjH/j) = ut- 
We denote the signal-to-noise ratio (SNR) p = P/a"^ and let R^; = PQ. The constraints in ([2]) now 
become Q ^ 0, Tr(Q) = 1. The secrecy rate maximization problem can be written as 

max _ CsiQ) = logdet(I„^ + pH^QH^j) - logdet(I„, + pH^QHy. (4) 
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According to Lemma [T] we assume hJjH/j — H^H e is positive semi-definite or indefinite (except the 
MISO channel). In llT4l . the authors assumed that hJ^H/j and H^jH^; are both positive definite, and 
hence they are both invertible. Here, we do not make that assumption. In fact, when > njj (and/or 
riT > tie), H^Hjj (and/or H^H^;) always have low rank and hence are not invertible. 
We also denote the feasible set of dJ]) as 

= {QIQ ^ 0,Tr(Q) = 1} (5) 

which is a convex set. 



III. Closed Form Secrecy Capacity of MISO Wiretap Channel 

We first provide a lemma that will be used here and in the proof of Theorem [6] later. 
Lemma 2: Let r and s be two known non-zero vectors, and rr^ — ss^ ^ 0. 

(i) If r = for a certain scalar ^, rr^ — ss^ has only one nonzero eigenvalue equal fo(|^p — l)||s|p 
with the associated eigenvector s/||s||; 

(ii) If r^s = 0, rr^ — ss^ has only two nonzero eigenvalues, i.e., rji = ||r|p, r]2 = — ||s|p with associated 
eigenvectors r/||r||, s/||s||, respectively. 

(iii) If neither r = nor r^s = 0, rr^ — ss^ has only two nonzero eigenvalues, i.e., rji = ||r|p — 
|c2||r^s| > 0,7^2 = ||r|p — |c4||r^s| < with the associated eigenvectors ei = c^^ ""^^^(r + 
|c2|e'*^'^~'^-'s), e2 = Cg ""^^^(rH- |c4|e''''^~''')s), respectively, where tp is the argument of t'^ s, i = \f—\. 



ci = ||r||2 + |c2|2||s||2 - 2|c2||rts|, |c2| = (||r||2 + ||s||2 - ^(||r||2 + ||s||2)2 _ 4|rts|2)/(2|rts|), 
C3 = ||r||2 + |c4|2||s||2 - 2|c4||rts|, |C4| = (||r||2 + ||s||2 + v/(||r||2 + ||s||2)2-4|rts|2)/(2|rts|). 
The proof is simple, therefore, omitted for the sake of brevity. But we outline the proof here. For the 
case (i), (ii), the proof is obvious. For the case (iii), first, we can show rr^ — ss^ has rank two, thus it has 
only two nonzero eigenvalues. Second, we assume the eigenvector has the form of a linear combination 
of r and s, and then show that this is indeed the case. 

Before discussing the general MIMO wiretap channel, we analyze a special case, i.e., the MISO wiretap 
channel in which the legitimate receiver and eavesdropper both have a single antenna, i.e., nji = ue = 1. 
Denote the channel vectors as hji and h^;. We give the following theorem. 

Theorem 1: The closed form expression for secrecy capacity of MISO wiretap channel is given by 

Cs = log (6) 

la 

where a = l + p\\h.E\?, b = 2 + /)||hij||2 + p||hs||2 + /32(||h^||2||h^||2 _ |htjh£;|2) and c = l + p||hij||2. 
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Proof: The secrecy rate maximization problem can be written as 



max C.(Q) = logi±4^ (7) 
Q>:0,Tr{Q)=l 1 + pWj^q,h.E 



which is a fractional program EOl associated with the following parametric problem 

F{a) = max [l + ph^oQh/j - ail + /ohj^Qhs)! (8) 

QhO,Tr(Q)=l L " \ 

where a > 0. Let a° be the unique root of F{a) = 0. According to |20|, the optimal Q corresponding 
to F{a°) also optimizes ([7]). Based on the fact that h^^Qh = Tr(Qhh^) for any vector h, we rewrite the 
optimization problem ([H) as 



Fia) = max 

QbO,Tr(Q)=l 



1 - a + pTr{Q(h^jhJj - ah^^h^ } 



(9) 



By eigen-decomposition h/jh^ — ah£;h^ = UaDaUa and letting Qq, = UaQUa, we obtain Qq, ^ 0, 
Tr(Qa) = 1, Q = U^Qc^uL It holds 

Tr{Q(hi^h]^ - ahehyi = TV{Q„D„} = TV{diag(Q,)D„} < X^U^Rh^R " "h^^h^ . (10) 

Equation (ITOl ) holds with equality if Qq is diagonal and has a unique nonzero entry (equal to one) 
corresponding to position of the largest entry in Dq. In other words, Q and h^jh^j — ah^jh^ have 
the same eigenvectors, and Q has rank one. Thus, it holds Q = maxUa,max where max is the 
eigenvector associated with the largest eigenvalue of hjjhjj — Qh£;h^. The largest eigenvalue and the 
associated eigenvector of h/jh|j — ah^jh^ can be expressed in closed form based on Lemma 2. In our 
problem, r = h.R, s = ^/ah.E■ By using Lemma |2l we now can obtain 

F(a) = 1 + ^^^- fl + ^^^^a + gV(l|hHP + a||hsP)2-4a|hihE|2. (11) 



. ^ , , ^ , \/lllllt^ll~-t-(J'llllfv'll~J~ H-fJ^II*"^ 

2 V 2^2 

F{a) = has a unique root given in closed form: 

a = . (12) 

2a 

The optimal Q is given by Q° = eie| where ei is defined in Lemma |2] where r = h^, s = \'7y°h.E- 
The secrecy capacity is given by Cs = loga°. 

■ 

Based on Theorem [B if = ^h^ and |^| < 1, then b = a + c, a — c>0 and further Cs = 0. 
This is consistent with the fact that when the legitimate receiver channel is a degraded version of the 
eavesdropper channel the secrecy capacity is zero. If h.R = ^Yie and |,^| > 1, then b = a + c, a — c<0 
and Cs = log((l + />|CPl|h£p)/(l + p||h£;p)) > 0. This is consistent with the fact that when the 
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eavesdropper channel is a degraded version of the legitimate receiver channel the secrecy capacity is 
positive. If hij ^ O^e, then b > a + c and it always holds that Cg > 0. Thus, if h.R / ^h^;, the MISO 
wiretap channel always has positive secrecy capacity independent of the channel. 

To gain more insight into the secrecy capacity, we consider the rate at which the secrecy capacity 
scales with logp as in |15|. If h/j / ^h^;, then under high SNR, it follows from ^ that 

Csip) = logp + logdlhnf - \h^j,hE\V\M^ + 0{l/p)) (13) 

where O(-) is the big-0 notation. The secrecy degree of freedom (s.d.o.f.) (also see (15)) of the MISO 
wiretap channel is given by 

s.d.o.f = lim = 1. (14) 

p^oo log p 

IV. Conditions for Optimal Input Covariance Matrix of MIMO Wiretap Channel 

In this section, we analyze a general MIMO wiretap channel. First, we obtain the necessary conditions 
for the optimal Q by using Karush-Kuhn-Tucker (KKT) conditions. Let us construct the cost function 

L(Q, d, *) = C,(Q) - 0(Tr(Q) - 1) + TY(*Q) (15) 

where 6 is the Lagrange multiplier associated with the constraint Tr(Q) = 1, 'I' is the Lagrange multiplier 
associated with the constraint Q h 0. The KKT conditions enable us to write ll2ll 

- ein^ + * = 0, (16) 

* ^ 0,Tr(*Q) = 0,Q ^ 0,Tr(Q) = 1, (17) 

where 

= pHtj(I„^ + pHKQH^)-^Hfi - pH^,(I„, + pUEQn^Ey^UE. (18) 



Here we use the facts: ^Tr(*Q) = and 



dQ 

For future use, we also rewrite ([T8] ) as 



aiogdet(I„„ +pH«Ql4) ^ [^hUi., +pH.QHt )-H.]- (19) 



= S^j(I„, + QSr)-^ - SE{In^ + QS^;)-^ (20) 

where Sr = /oHJ^H^j, = pH^Hs, which follows from the fact: H(I + pQHtH) = (I + /)HQHt)H 
for any matrix H and hence (I+pHQHt) = H(I + pQHtH) i, Ht(I+pHQHt) = HtH(I + 
pQHtH)-!. 
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In this paper, is an important variable for the optimal input covariance problem. It has the following 
property which will used later. 

Property 1: For any Q h 0, Amax(0) > 0; For any Q ^ 0, Tr(Q0) < Tr(Q)Ai„ax(0), and in 
particular, for any Q G fi, Tr(Q0) < Aniax(0)- 

The proof is given in Appendix |B] 

From the KKT conditions ([T6l l and ([TT] ). we obtain the equivalent (but without containing the Lagrange 
multipliers) conditions for optimal Q consisting of a set of equations given in the following theorem. 
Theorem 2: The optimal Q ^ satisfies 

Q0 = Tr(Q0)Q (21) 
Amax(e) = Tr(Q0). (22) 

Please see Appendix ICl for details. 

Equations (|2TI ) and (l22l ) provide two elementary conditions that characterize the optimal Q. At this 
point we do not have a proof that any Q satisfying the conditions of Theorem |2] is the optimal input 
covariance. However, for some special cases, e.g., the MISO wiretap channel analyzed in ^IIll this is 
true. In particular, for this case we provide the following theorem. 

Theorem 3: For MISO wiretap channel, any Q satisfying the conditions of Theorem^ is the optimal 
input covariance. 

The proof is given in Appendix |D] 

Now we proceed. From Property [U and (l22l ). we know that the optimal Q satisfies 

TY(Q0) > 0. (23) 

Based on (|2TI ) and (1231 ). and by taking trace operation over both side of (|2TI ). it can be easily seen that 
Tr(Q) = 1. That is to say, equations (EB and ^ imply Tr(Q0) > and Tr(Q) = 1. 

The condition (|2TI ) reveals that the optimal Q satisfies that Q and commute and have the same 
eigenvectors ll22l p.239]. The condition (l22l ) means that the eigenvalues of corresponding to the positive 
eigenvalues of Q are all equal to Tr(Q0), while the remaining eigenvalues of (i.e., corresponds to 
the zero eigenvalues of Q) are all less than or equal to Tr(Q0). Obviously, if the optimal Q has full 
rank, then = for a certain > 0. 

It can be shown that based on the conditions of Theorem |2j the optimal Q has the following properties. 

Property 2: The optimal Q satisfies: 
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(i) rank(Q) = rank(Si? — S^;) > rank(Q); 

(ii) Q(Sr - Si?)Q ^ 0; 

(iii) Q + QS£;Q and Q + QS/jQ commute and have the same eigenvectors. 

For readability, we put the proof of Property |2] in Appendix El A direct result of Property |2] is the 
following: 

Property 3: when Sr — Se ^, the optimal Q has low rank. 
The proof is simple. When S/j — S^; is indefinite, if the optimal Q has full rank, then Property |2] (ii) 
leads to — S^; ^ 0, which violates that Sr — S^; is indefinite. When Sr — S^; ^ but S/j — S^; ^ 0, 
it follows from Property |2] (i) that the optimal Q has low rank. This result was also pointed out in iTHI . 
When Sr — Se ^, the optimal Q may have low rank or full rank. 

Before ending this section, we point out that we can combine the elementary conditions (|2TI ) and (l22l ) 
into a single equation. When © and Q commute and have the same eigenvectors, + ^Iut ^i^^ Q 
commute and have the same eigenvectors for any real number 7, and vice versa. We can find a certain 
7 such that + 7lnr >~ for any Q ^ 0. Based on in (l20l) . we have 



Amm(0) > -A max (Se). (24) 
When 7 > Aniax(S£;), it always holds that + ^\nT 

^ 0. Let K = + 7!^^ and hence K ;^ 0, 
Tr(QK) > for any Q ^ but Q / 0. Equations dH) and (l22]l are equivalent to 

QK = Tr(QK)Q, (25) 

An,ax(K) = Tr(QK). (26) 
We can combine the above two equations to a single one as follows. 

QK = i (^TV(Q) + K.U^n (27) 

or equivalently, 

QK=i(TY(Q) + ^)Q (28, 

where K = K/Amax(K). We give the following theorem. 

Theorem 4: Any Q ^ that satisfies H27[ also satisfies the conditions of Theorem^ 

Proof: It is easy to verify that i(Tr(Q) + l/Tr(Q))Amax(K) is an eigenvalue of K. But Amax(K) 
is the largest eigenvalue of K. On the other hand, i(Tr(Q) + l/Tr(Q)) > 1 holds with equality if and 
only if Tr(Q) = 1. Thus, we know Tr(Q) = 1. With this, taking trace operation over both sides of ( |27l ) 
leads to Ai„ax(K) = Tr(QK). ■ 
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Summarily, we can alternatively do one of the following two things to find Q that satisfies the conditions 
of Theorem ID 

(i) Find Q ^ satisfies ^ and (Egjl; 

(ii) Find Q ^ satisfies (|27] ): 

We will discuss the algorithm to search for such Q in Will 

In the following sections, we will analyze some special cases. In particular, for n-r = 2 we obtain the 
optimal Q in closed form. If H^H^ — H^H^; has all negative eigenvalues except a positive eigenvalue, 
we show that the optimal Q has rank one and can also can be expressed in a closed form. For general 
cases, we prove that the optimal Q is a fixed point of a mapping from a convex set to itself, and propose 
an algorithm to search for it. 

V. The Case ut = 2 

In this section, we analyze the case ut = 2, i.e., the transmitter has two antennas. It includes the four 
cases {riT, nR.UE) = (2, 2, 2), (2, 2, 1), (2, 1, 2), (2, 1, 1). In gml the MISO wiretap channel with nx = 2 
belongs to (n^, nR,nE) = (2, 1, 1). In [ 13], the case {tit, nji,nE) = (2, 2, 1) is analyzed. We derive the 
optimal Q in two subsections in which S/j — S^; ^ or Sr — S^; ;^ 0. We also analyze the rank of 
optimal Q with respect to SNR. 

A. Sr-Se^O 

According to Property [3l the optimal Q has low rank (rank one) and hence it has the form Q = uu^ 

where u is a unit-norm vector to be determined. We can rewrite 

n fn,\ 1 l + u1'S/ju ut(l2 + Sif)u 

Q = log , — = log ^„ . • (29) 

l + uTS£;U uT(l2 + S£;)u 

The optimal Q is easily obtained to be Q° = u°u°^ where u° is the eigenvector associated with the 
largest eigenvalue of (I2 + S_e)^^(I2 + Sr). The secrecy capacity is given by 

Cs = log (A^ax{(l2 + SErHh + Sr)}) . (30) 

We can express Cg in closed form. Denote 

H^H£;= I °i ,HtjH^= I ' " I . (31) 
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By using the fact: for any 2x2 matrix A with two real eigenvalues, the largest eigenvalue is given by 
Amax(A) = [Tr(A) + V(Tr(A))2 -4det(A)]/2, and the matrix inverse formula 



an ai2 \ 1 



ai2 a22 J 011022 -|ai2P , -a^^ 




(32) 



we can obtain 



r - 1 A + ^A'^-4{B,-B2) 

- 2[(l/p + ai)(l/p + ci)-|6iP] ^^^^ 

where A = aiC2 + a2Ci -&162 -^1^2 + (ai + a2 + ci + C2)/p + 2/^^, Bi = (a2Ci -6162 + (02 + ci)/p + 

l/p^)(aiC2 - blb2 + (ai + C2)/p + and B2 = (62C1 - 61 C2 + (62 - h)/p){aib* - 026^ + {b* - b\)/p). 

Now we analyze the secrecy degree of freedom (s.d.o.f.) defined in (fT4l ) which is different whether 

S E has full rank or not. 

• Case 1) S^; has rank two (full rank) 

In this case, noting that (I2 + Si5)-i(l2 + S^) = (h/p + U^EUE^^ih/ P + HrHr), we have 
Cs log(Amax{(H^H£;)"^HjjH/j,}) as p ^ oo. Thus, we get 

s.d.o.f = lim = 0. (34) 

p^oo log/9 

• Case 2) S^; has rank one (low rank) 

In this case, H^H^; is singular, hence can be expressed as H^H^; = V2V2. By using (l60l ). we can 
write (I2 + pS2)-ni2 + pSi) = pih - pv2vt/(l + p||v2f ))(l2/p + Si) ^ pih - S2/Tr(S2))Si 
as p ^ 00. Thus, as p ^ 00, Xmi,x{{h + pS2T^{h + pSi)) /)Amax((l2 - S2/Tr(S2))Si). Thus, 
we get 

s.d.o.f = lim = 1. (35) 

p-»oo log p 

We can also use (l33l) to obtain the same result. 
B. Sr-Se^O 

In this case, the optimal Q may have fuU rank or low rank. If the optimal Q has low rank, it is given 
in ( [29l ). Therefore, in the following we focus on the case in which the optimal Q has full rank. The 
optimal Q can be determined from the above two cases. 

Since Q ^ 0, it follows from (|2TI ) that G must be a positive scalar multiplication of I2. Recall from 
(I20I ) that e = Snih + QSny^ -SEih + QSE^^. We know Sr >- 0, but S^; is not necessarily positive 
definite. Thus, in the following, we discuss two cases respectively: a) S^; has rank two (full rank); b) 
S^; has rank one (low rank). 
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• Case a) S^; has rank two (full rank) In this case, ;^ 0, S^; ;^ 0. We can rewrite 

= (S^i + Q)-i - (S^i + Q)-i = 9h,9 > 0. (36) 

Based on the eigen-decomposition (S^^ — S^^)/2 = UiDiUj where Di = diag((ii, ^2), di > 
d2 > but d\ + d\ ^ (otherwise, Sr = violates the assumption 7^ S^;), and letting 
C = UJ(S^^ + S^^)Ui/2, we get S^^ = Ui(C + Di)u| and S^^ = Ui(C - Di)ul. On 
inserting the latter expressions in 061 ) we get 

(Q + C - Di)-i - (Q + C + Di)-i = ^ > (37) 

where Q = u|QUi. Note that Di is diagonal. We can actually show that Q + C must be diagonal. 
To prove this, let us denote the (1, 2)th entry of Q + C by 512- We know that the (1, 2)th entry of 
(Q + C - Di)-i - (Q + C + Di)-i equals 

- gi2[det(Q + C - Di)-i - det(Q + C + Di)"^] = (38) 

which leads to qi2 = 0. Here we used ([321) and the fact det(A) > det(B) for A ^ B, B ^ 0. 
Since Q + C is diagonal, we denote Q + C = diag(gi,g2) and Q = diag(gi,g2) — C. It follows 
from Q ^ 0, Tr(Q) = 1 that + 52 = 1 + Tr(C), diag(gi, ^2) >- C. Combining these with ([371) 
results in 

(diM - l)ql + 2(1 + Tr(C))g2 + di{di - ^2) - (1 + TY(C))2 = 0. (39) 

We can solve 52 from the quadratic equation ( [39l ) and qi = \ + Tr(C) — q2- If diag(qi,g2) >- C 
holds, then Q = UiQUj is a possible solution. If the equation ( [39l ) has no positive roots or 
diag( 51,^2) ^ C, it means the optimal Q has low rank. 

• Case b) S^; has rank one (low rank) In this case S^; can be expressed as = Vgvl. We eigen- 
decompose S^; = Uediag(Ae, 0)Ue. Similarly, we get 

= (S^i + Q)"i - Sij(l2 + QSe)-i = ei2, ^ > 0. (40) 

Let us define Q = uiQUe, Sr = UeS^^Ue. Inserting Q and Sr into ( |40l ) results in 

^1 Ap / 1 \ 

(Sii + Q)-^-— ^ \=ei2,e>Q. (41) 

1 + qiiK \ y 

where is the (1, l)th entry of Q. It follows from (|4TI) that Sr + Q is diagonal. Thus, we denote 
Sfl + Q = diag(gi,g2) and Q = diag(gi,g2) - ^R- It follows from Q ^ 0, Tr(Q) = 1 that 
9i + 92 = 1 + Tr(Sij) and diag(gi, ^2) >- S/j. Combining the above and (|4TI) results in 

Kql + 2(1 - (SH,)iiAe)gi - (1 - (Sfi)nAe)(l + 1^(8^)) = (42) 
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where (Si?)ii is the (1, l)th entry of Sr. We can solve qi from the quadratic equation (l42l ). and 
then get q2 = 1 + Tr(S/j) — qi. If diag(gi, 52) >- S/j holds, then Q = UeQUl is a possible solution. 
If the equation (|42)) has no positive roots or diag(gi,g2) ^ Sr, it means the optimal Q has low 
rank. 

C. Rank of Optimal Q 

For the non-wiretap MIMO channel the rank of optimal input covariance has a non-decreasing property 
with respect to SNR |[T9]| . In this section we consider the behavior of the rank of optimal input covariance 
of the MIMO wiretap channel with respect to SNR. 

When Sr — Se ^ 0, according to the result in W-A[ the optimal Q has rank one, independent 
of SNR, and hence follows the non-decreasing property of rank. Next we focus on S/j — S^; >- 0. 
According to W-A[ if the optimal Q has rank one, it can be expressed as Q = uqUq where uq is the 
eigenvector associated with the largest eigenvalue Aq of (I2 + pS2)~^(l2 + pSi) where Si = hJ^Hr, 
S2 = H^Hs. Denote Sq = Si(l2 + puouJSi)"^ - 82(12 + puouJS2)"^. Since Si - S2 ^ 0, we can 
rewrite Sq = (I2 + pSiUoUq)~^(Si — S2)(l2 + /OUoUqS2)^^ and state that Sq has full rank (rank two). 
Denote 

^ ^ ^ t« rj. 4(2Si + pSf)uo 4(2S2 + pSi)uo 
g{p) = 2uiSoUo - Tr(So) = ^p- — j- — Tr(Si - S2) (43) 

+ pbi jUo U|!)(i„^ + /9b2jUo 

where the matrix inverse formula ( [6OI ) is used. We give the following result. 

Lemma 3: If g{p) < 0, then the optimal Q has rank two; If the optimal Q has rank one, there must 
be g{p) > 0. 

Please see Appendix |F] for details. 

Then we can prove the following result. 

Theorem 5: limp_^oo g{p) < 0, hence, according to Lemma\3\ there exists a certain pQ such that when 
p > pq, the optimal Q has rank two. 
The proof is given in Appendix |Gl 

Theorem [5] reveals that when the SNR is sufficient large, the optimal Q always has rank two. At this 
point, we do not prove the rank non-decreasing property of the optimal Q for the case Sr — Se >- 
completely. 

VI. Sr - Se Has All Negative Eigenvalues Except One Positive Eigenvalue 

We analyze the case in which hJjH/j — H^H^; has all negative eigenvalues except one positive 
eigenvalue, e.g., ut = 3, H^H/j — H^H^; has two negative eigenvalues and a positive eigenvalue. In 
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particular, this always occurs when Se has full rank and nR = 1, i.e., the legitimate receiver has a single 
antenna, as the following lemma stated. 

Lemma 4: Let x is a known non-zero vector, X is a known positive semi-definite matrix. Assume 
xx^ — X is indefinite or positive semi-definite. If X hasfiill rank, or X has rank nx — 1 and x is linearly 
independent of the eigenvectors associated with the non-zero eigenvalues of X, then xx^^ — X has all 
negative eigenvalues except one positive eigenvalue. 

Proof: First, we prove the case that X has full rank. Let A be any eigenvalue of xx^^ — X. It holds 
det(AI - xxt + X) = 0. When A > 0, noting that AI + X ^ 0, we get det(AI - xxt + X) = det(AI + 
X) det(I-(AI + X)-ixxt) = det(AI + X)(l-xt(AI + X)-ix) = which leads to xt(AI + X)-ix = L 
Here we use the fact det(I + AB) = det(I + BA). It is easy to prove that x^(AI + X)^^x decreases 
strictly with A. Thus, there is at most one A > such that x1^(AI + X)^^x = 1. This, when combined 
with the fact that xx^ — X is indefinite or positive semi-definite, gives the desired result. 

Second, we prove the case that X has rank n-p — 1 and x is linearly independent of the eigenvectors 
associated with the non-zero eigenvalues of X. Denote the eigen-decomposition X = Aiuiu| + • • • + 
\nT-i^nT-i'^\r^-v We can write xx"!" - X = F^D^fJ. where F^; = [x, ui,--- ,u„^_i], D^; = 
diag(l, — Ai, • • • , — A„y_i). Since x is linearly independent of the eigenvectors associated with the non- 
zero eigenvalues of X, it holds that F^; has full rank. According to Sylvester's law of inertia 1231 p.223], 
we know F^^D^^F^. and D^^ have the same number of positive, negative, and zero eigenvalues, thus the 
desired result is obtained. 

■ 

According to the above lemma: when n/j = 1 (hence, Sr can be expressed = v^vj) and S^; has 
full rank, ^r — Se has all negative eigenvalues except one positive eigenvalue; when n/j = 1, S^; has rank 
nj' — 1 and is linearly independent of the eigenvectors associated with the non-zero eigenvalues of S e, 
Sr — Se has all negative eigenvalues except one positive eigenvalue. But we point out that it does not 
limit to the cases in Lemma |4] in which hJjH/j — H^H^; has all negative eigenvalues except one positive 
eigenvalue. In fact, this will even occur when Sr and S^; both have full rank. When H^H/j — H^H^; 
has all negative eigenvalues except one positive eigenvalue, we give the following theorem. 

Theorem 6: If H^H/j — H^H^; has all negative eigenvalues except one positive eigenvalue, the 
optimal Q has rank one. Also, the optimal Q is given by Q° = u°u°^ where u° is the eigenvector 
associated with the largest eigenvalue of (I^^ + S£;)^^(I„^ + Sr). The secrecy capacity is given by 

Cs = log (A„,ax{(In. + SE)-Hlnr + Sr)}) . (44) 
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Please see Appendix iHl for details 

4 



Similar to if HtHi? has full rank, we obtain 



s.d.o.f = lim = (45) 

and if H^Hg has low rank, we obtain 

s.d.o.f = lim = 1. (46) 

p^oa log p 

VII. Algorithm for General MIMO Wiretap Channel 

In this section, we propose an algorithm to search for the optimal Q which applies for any MIMO 
wiretap channel. The algorithm is based on the conditions of Theorem |2] or Theorem |4l 
It follows from (US]) that 

k1/2qkV2 ^ 'iv(QK)Q (47) 

which enables us to get 

K1/2qk1/2 

'^"^yir*^"^'^ '''' 

Note that /(Q) ^ and Tr(/(Q)) = 1 for any Q G fi. The equation (l48l) defines a mapping from a 
convex set to itself: — > 17, Q /(Q)- The optimal Q corresponds to a fixed point of /(Q), i.e., 
/(Q°) = Q°- To search for the fixed point, the iterative expression is 

Q'+' = /(Q'),fe = 0,l,--- (49) 

The initial point Q° can be set to I„y or choose a good initial point. The iterations stop when HQ'^"''^ — 
Q*^|| < 10^®. If the convergent Q satisfies ( [261 ). we obtain a solution satisfying the conditions of Theorem 
[21 otherwise, we choose a different initial point. 

VIII. Numerical Simulations 

We give some examples to illustrate the proposed algorithm. For illustration purpose, we consider a 
MIMO wiretap channel where = 4, n/j = 4, n^; = 3. 

First, we take an example for H^H/j — H^H^; >- 0. The channel matrices are given by 

/ -0.1107 -0.1225i 0.0582 - 0.3483i 0.3239 - 0.0071i -0.2872 - 0.2655i ^ 
0.5128 - 0.3239? -0.8903 - 0.0318i -0.5524 - 0.0365z -0.2072 + 0.3047i 

(50) 

-0.0041 + 0.0265i 0.0871 - 0.0253i 0.0183 + 1.1679i -0.0784 + 0.0415i 
y -0.4699 - 0.1014Z -0.0888 + 0.1127i 0.2099 + 0.3282i 0.1734 - 0.4146i y 



H 



R 
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and 



H, 



/ -0.0766 + 0.1370i -0.0977 - 0.0985i 0.0002 - 0.0695i 0.0583 + 0.0356i 
-0.0355 - 0.1167i 0.1607 - 0.109H -0.0809 + 0.1481i -0.0218 + 0.1109i 
Y 0.1375 - 0.0381i -0.0845 - 0.0610i -0.0011 + 0.1129i -0.0393 + 0.1124? 



(51) 



The eigenvalues of H^jHi? - H^H^; are 0.0085, 0.3704, 0.8945, 2.5213. Fig. [I]|3] depict respectively the 
eigenvalues of Q'^, secrecy rate and ||Q'^^^ — Q'^H in the iterations where the SNR is p = 8dB. Fig. 
|4]|5] depict respectively the (possible) secrecy capacity and eigenvalues of (possible) optimal Q under 
different SNRs. It can be seen from Fig. |5]that the (possible) optimal Q can have rank one to four with 
the increasing SNR, which shows that when h|jH/^ — H^H^; >- 0, the (possible) optimal Q may have 
low rank or full rank. 

Secondly, we take an example for hJjHjj — H^H^; ^ 0. The channel matrices are given by 

\ 



H 



R 



( -0.1110 - 0.0667i 
0.2877 + 0.6779i 
0.3266 - 0.2779? 



-0.1937 - 0.1349i -0.0752 - 0.2707i 
-0.7832 - 0.2249i 0.4350 + 0.2637i 
-0.2345 - 0.4472i 0.2448 + 0.3488i 



-0.2718 + 0.2730i 
0.4160 + 0.5109i 
-0.6794 - 0.0117i 



(52) 



y -0.1221 + 0.4915? 0.0959 - 0.2557i -0.0219 + 0.5077i 0.1449 + 0.3294? j 



and 



H, 



/ 0.1468 - 0.1185i 0.4071 + 0.4469? 

-0.0520 + 0.2917i -0.4978 + 0.0545i 
y 0.5799 - 0.1767? 0.2298 + 0.3331? 



0.2474 - 0.3291i 
0.0779 - 0.3472? 



-0.6264 - 0.1313? 
-0.0132 - 0.1327? 



-0.1151 - 0.2000? 0.1404 - 0.3501? 



(53) 



The eigenvalues of H^H^j - H^Hij are -0.8206, -0.1565,0.9365, 1.8506. Figures HM] depict respec- 
tively the eigenvalues of Q'^, secrecy rate and ||Q^^^ — Q^|| in the iterations where the SNR is p = 8dB. 
Figures l9lfT0l depict respectively the (possible) secrecy capacity and eigenvalues of (possible) optimal Q 
under different SNRs. It can be seen from Fig. [TO] that the (possible) optimal Q always has rank two, 
which equals the number of positive eigenvalues of h|jH/j — H^H^;. 



IX. Conclusion 

We have investigated the problem of finding the optimal input covariance matrix that achieves secrecy 
capacity subject to a power constraint. In particular, for the multiple-input single-output (MISO) wiretap 
channel, the optimal input covariance matrix is obtained in closed form. For general cases, we derive the 
necessary conditions for the optimal solution consisting of a set of equations. For the case in which the 
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transmitter has two antennas, the derived necessary conditions can result in a closed form solution. If the 
difference is indefinite and has all negative eigenvalues except one positive eigenvalue, we prove that the 
optimal input covariance matrix has rank one and can be obtained in closed form. For other cases, we 
prove that the solution is a fixed point of a mapping from a convex set to itself and provide an iterative 
procedure to search for it. 

Appendix A 
Proof of Lemma [T] 

First we show the necessary part. When H^jH/j— H^H^; is negative semi-definite, so is Ry^(H^H^— 
H^H£;)Ry^ which leads to In^ +kI^^'H.^j^'H.e'RI^^ /a^ >z +Ry^H^HRRyV^^- With this, using 
the fact: if A ^ 0, B ^ 0, A ^ B then det(A) > det(B), and applying the identity det(I + AB) = 
det(I + BA) to ^ results in Cs(R^) < 0. To show the sufficient part, we rewrite ([3]) as 

Cs(Rx) = logdet (in^ + ^Ry2(HtjHR - li^j^liE)Ii'J\ln^ + Ry^H^HeRyV^^")"') • (54) 

Note that + Hx^'^H^^JIe^x^'^ / cr^)^^ >- 0, and it suffices to show that there exists R^. such that 
Ry^(HjjH/j — H^H£;)Ry^ >: 0. Let us define the eigen-decomposition H^jH^j — H^H^; = UrDrUj. 
It is easy to verify that R^^ = UrD^^uJ is a choice where the entries of the diagonal D^; are zero 
corresponding to the position of negative entries in D^. 

Appendix B 
Proof of Property [T] 

First, we prove the former part. If Sr >- and S^; ;^ 0, then we can rewrite ( [201 ) as 

e = (S^i + Q)-i-(S^i + Q)-^ (55) 

We can state that is not negative semi-definite, otherwise, we get S^; ^ Sr which violates the 
assumption that Sr — S^; is indefinite or positive semi-definite. 

Next we consider the case that Sr or S^; are singular. Denote Sr — S^; = A. By using the fact: 
Se{1 + QSe)~^ = (I + SeQ^^Se, we can rewrite ^ as 

= SR{In^ + QSr)-^ - (I„, + SeQ)-^Se 
= (I„, + SeQ)-'A{I^, + QSr)-^ 

= {In^ + SijQ)-ip(In, + QSe)-^ (56) 
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where P = A(I„^ + + QSb)"^QA)"^ According to Sylvester's law of inertia 1231 p.223], it 
suffices to prove that P has positive eigenvalue. If A is nonsingular, we can write 

P = (A-i + (I„, + QSi,)-iQ)-i. (57) 

Note that (I„^ + QSe)'^Q >z and using the assumption A is indefinite or positive semi-definite, 
hence, we get that P has positive eigenvalue. If A is singular, there exists > such that for any 
0<(5<(5o, A' = A — 6InT is nonsingular and also indefinite or positive semi-definite. Similarly, we 
can prove that P' = A'(I„y + (I^^ + QS£;)^^QA')^^ has positive eigenvalue. Next we prove P >- P'. 
Denote W = (I„^ + QS£;)~^Q. Similar to the skill in ([561), we get 

P-P' = (I„, +AW)-iWi(I,,+WA)-i (58) 

where Wi = 6{Im, —6{In.j, + WA)~^W)^^. We can see that when S is sufficient small, Wi is positive 
definite. Thus, P ;^ P' and P has positive eigenvalue. 

Second, we prove the latter part. Denote the eigen-decomposition = UeDeUg) and let Qi = 
U^QUe. We know Qi ^ and Tr(Qi) = Tr(Q). With these, we can write 

Tr(Q0) = Tr(QiDe) = Tr(diag(Qi)D0) < Tr{Qi)X^,,(De) = Tr(Q)A,nax(0). (59) 
In particular, if Tr(Q) = 1, we get Tr(Q0) < Amax(©)- 

Appendix C 
Proof of Theorem [2] 

It follows from ([TT] ) that ^'Q = Q^' = 0, that is, ^' and Q commute and have the same eigenvectors 
II22I p.239] and their eigenvalue patterns are complementary in the sense that if Aj(Q) > 0, then Xi{^) = 
0, and vice versa {Tf]. This result, when combined with (fT6l) . implies that and Q commute and have 
the same eigenvectors, i.e., they have the eigen-decompositions Q = U^DgUj and = UqDeUj. 
Further, we get 0Q = Q0 = ^Q, which, when combined with Tr(Q) = 1 and the fact Tr(Q0) = 
Tr(Qi/20Qi/2) is always real, leads to 6" = Tr(Q0) and (EU) (also see [18]). 

The condition (|2TI ) reveals that for the optimal Q, Q0 is a scaled version of Q. Further the eigenvalues 
of corresponding to the positive eigenvalues of Q are all equal to Tr(Q0), while the remaining 
eigenvalues of are all less than or equal to Tr(Q0), which follows from (fT6l ). ((2T]) and ^ 0. Based 
on the above, it holds the second condition (l22l ). 
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Appendix D 
Proof of Theorem [3] 

In the MISO wiretap channel, S/j = ph/jhj^ and S^; = ph£;h^. By using the matrix inverse formula 
for two vectors x and y 

(I + xyt)-l =I-xyt/(l + ytx), (60) 



We can write 

1 + pWj^QhR 1 + ph^Qhe ■ 
That is to say, has the form of aih/jhjj — a2h£;h^, ai > 0, 02 > 0. According to Lemma [2l we 
know: if = ^h^;, then has only one nonzero eigenvalue; if h/j / h^;, then has only two nonzero 
eigenvalues, one is positive and the other is negative. With this, since Q satisfies (|2TI) . it is easy to verify 
Q has rank one. Let Q = uu^ and we have 

^ R R H E E (^2) 

1 + ph^uuth^ 1 + ph^uuth^; 

IV(Q0) = 1 1 . (63) 

1 + ph^uuth^; 1 + ph|juuthj:j 

Let wi = 1 + ph^uu^hij, u;2 = 1 + ph^uu^h^;. According to Lemma |2l the largest eigenvalue of is 
given by 



ZWl ZLJ2 ^ V V "^1 (-^2 / OJ1UJ2 



Since Q satisfies (I22b . we have 



p\\hR\\^ pWhE^ ^ 1 ( P\M' ^ p\\hE\\y _ V|hl;h^P ^ J_ _ J_ (65) 



which leads to 



_/',_,_ P\\^Ef\ ^^P, { „. ,12 ^ ||u„||2'ilV _ .lU 



l + ::^_,l + ^,_^ + ^WMM. + ||hdP-.l _4|hi,h,p-^=0. (66) 



This equation ( |66l l is exactly (fTTl) . i.e., F{a) = where a = uJxjbj^- On the other hand, we know 



wi 1 + phjjuu'l'hij 
^2 1 + phtuuth^; 



(67) 



According to the result in ^IIll the root of F(a) = corresponds to the maximization of the right 
hand side (RHS) of (|67] ) (see also (|7])). Thus, the conditions of Theorem |2] guarantee the optimal input 
covariance. 
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Appendix E 
Proof of Property [2] 

From dinil, we know that Q0 = (I„y + QS^)"^ - (I„y + QSr)"^. With this, Left-multipUcation by 
IriT + QSe and right-multipUcation by + QS/^JQ of both sides of (|2T] ) results in 

Q(Sfi - Se)Q = TY(Q0)(Q + QSijQ)(Q + QSrQ). (68) 

Note that the left hand side (LHS) of (l68]l is Hermitian, hence, (Q + QS£;Q)(Q + QS/jQ) is Hermitian 
as well, which implies that the matrices Q + QS£;Q and Q + QS/jQ commute and have the same 
eigenvectors. On the other hand, the RHS of ( [68l ) has all non-negative eigenvalues, thus Q(S/j — S£;)Q ^ 
0. 

Finally, recalling @ = (I„^+S£Q)"^(Sj?-S£;)(InT+QSR)"\ we know rank(0) = rank(SR- 
S^;) which follows from the fact rank(A) +rank(B) — n< rank(AB) < min{rank(A), rank(B)} for 
two nx n matrices A and B p3l. Further, since rank(Q0) < min{rank(Q), rank(0)}, if rank(0) < 
rank(Q), then rank(Q0) < rank(Q), but this violates (|2T] ). Thus, rank(0) > rank(Q). 

Appendix F 
Proof of Lemma [3] 

First, we can prove that ug is the eigenvector of Sq- To prove this, it suffices to show that SqUq = ?/oUo 

for a certain scarlar r]Q. By using the formula ( [60l ). we can write 

Siuo S2U0 
Squo = 7 1 . (69) 

1 + pU|^SiUo 1 + /)U^S2Uo 

Since uq is the eigenvector of (I2 + pS2)~^{l2 + pSi) associated with the eigenvalue Ao, it holds 
(I2 + pSi)uo = Ao(l2 + pS2)uo which leads to two facts: 1 + pugSiUo = Ao(l + /OUqS2Uo) and (Si — 
AoS2)uo = ^^'^uq. Inserting them into ( l69l ) gives SqUq = i]oUo where rjo = (Aq — l)/(/o(l + /ouJSiUo)). 
We also know rjo = UqSoUq. 

Second, let v is any unit-norm vector and v / ^uq. Define the function 

v) = logdet(l2 + p{{l - t)uouJ + tW)Si) 

-logdet(l2 + p((l-t)uouJ + ivvt)S2), tG[0,l]. (70) 

It is easy to know: (1 — t)uoUQ + tvv'l" G il; when t = 0, {1 — t)uoUQ + tvv"!" = uquJ; when t > 0, 
(1 — t)uouJ + tvv^ has rank two. We can state that if §f |t=o > for a certain vq, the optimal Q has rank 
two. The reason is simple: assume that the optimal Q has rank one, then it holds g{0, v) > g{t, v), Vt, v. 
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But on the other hand, it follows from §f |t=o > that there exists to > such that fi'(to) vq) > ^(O, vq). 
This produces a contradiction. By using the derivative formula |[T6l 

^»A±^=l..{B(A + <B)-'}, (71) 

we can obtain 

^lt=o = Tr{p(vv^ - uouJ)So} = p(v^Sov - uJSqUo). (72) 

It follows from (1721 ) that if UqSqUo < Amax(So), there always exists v such that §f|t=o > 0, and 
hence the optimal Q has rank two; if the optimal Q has rank one, there must be ^\t=o < which 
leads to UqSqUo = Amax(So)- Further, since uq is the eigenvector of Sq, thus UqSqUo < Amax(So) and 
uJSqUo = Amax(So) are equivalent to 2uqSoUo < Tr(So) and 2u|)SoUo > Tr(So), respectively. This 
completes the proof. 

Appendix G 
Proof of Theorem [5] 

We discuss two cases respectively. 

(i) S2 has full rank. Since (I2 + pS2)"^(l2 + pSi) = {h/p + S2)"Hl2/p + Si), we know that as 
p — > 00, uo goes to the eigenvector uq associated with the largest eigenvalue A of 82^^81. We can 
write Siuo = A82U0, hence UQ8fuo = AuqSiS2Uo, Uq8iUo = AuqS2Uo. Then 

uj8fuo U082U0 uJ(Si - 82)8220 rp /o c ^ 

g{p) g ~ - ~ Q ~ - Tr(8i - 82) = , ^ Tr(8i - 82 

U081U0 U082U0 U582U0 

< A^ax(Si - 82) - Tr(8i - 82) < (73) 

where we use the fact x^Ax/(x'''Bx) < Amax(AB^^). 

(ii) 82 has rank one and hence can be expressed as 82 = V2V2. As p ^ 00, 

uj8fuo UJ82U0 , uj8^uo rp /o \ rp /c c ^ 

9{p) -J- f- Tr(8i - 82) = -f- Tr(82) - Tr(8i - 82) 

u^biUQ U|!)b2Uo u^biuo 

< A„,ax(Si) - Tr(8i) < 0. (74) 

Appendix H 
Proof of Theorem [6] 

According to Property |2l we know Q(8ij — 8e)Q ^ 0. We will show that under the condition of 
Theorem |6l any Q ^ that satisfies Q(8/j — 8£;)Q ^ has rank one. To prove this, we use the 
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mathematical induction which consists of two steps: (1) showing that the statement holds when = 2; 
(2) showing that if the statement holds for some nr = k > 2, then the statement also holds when k + I 
is substituted for k. 

First, we address (1). It follows from Property [3] that the optimal Q has low rank, i.e., rank one. Next 
we deal with (2). We denote the eigen-decomposition Sr — Se = UDU^ where D is a diagonal matrix 
with all negative diagonal entries except a positive one. Let X = U^^QU. Then Q(S/j — S£;)Q ^ is 
equivalent to XDX >z 0. It suffices to prove the following: 

Problem: Assume the following is true: any k x k diagonal matrix D which has k — I negative 
entries and one positive one in its diagonal, then any kxk matrix X ^ that satisfies XDX >z 
always has rank one. Is it true for k + 17 
Denote the {k + 1) x {k + 1) diagonal matrix 



where Di is a A; x A; diagonal matrix with all negative diagonal entries except a positive one, ^ > 0. We 
also denote the {k + I) x {k + 1) positive semi-definite matrix 



where Xi is a k x k matrix, b is a x 1 vector, x is a scalar. It follows from X ^ that: Xi ^ 0, 
a; > 0; if a; = 0, then b = 1231 . Now we can write 



When X = 0, there will be b = and hence XDX = diag(XiDiXi, 0) ^ which is equivalent to 
XiDiXi >: 0. With this, based on the assumption for k, Xi has rank one, thus, X has rank one. In 
the following, we consider x > 0. Combining dTTl ) with XDX >z results in XiDiXi — ^bb^^ >z 
and b^Dib - > El. From b+Dib - ^x^ > 0, we get b^Dib > > (hence b / 0). 
It follows from XiDiXi - ^bbt ^ that XiDiXi >z 0. With this, by using the assumption for k, 
we know Xi has rank one and can be expressed as Xi = vv^^. Thus, XiDiXi = (v^Div)vv^ and 
v^^Div > 0. We can prove that v = cib for a certain scalar ci. To see why this is the case, let us rewrite 
XiDiXi-^bbt = (vtDiv)vvt-^bb^ Note that v^^Div > 0, hence according to Lemma|2l if v / cib, 
then (vtDiv)vvt — ^bb^ always has a negative eigenvalue. But this violates XiDiXi — ^bb^ >z 0. 




(75) 




(76) 




(77) 
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Thus, V = cib. Next, we determine ci. Now we can write XiDiXi — ^bb^ = (bTDib|ci|^ — ^)bb1' 
and 

/ (/3|ci|4-^)bbt (/?|ci|2-ex)b \ 
XDX = ^'"1 I ''^ vA-i I s / ^^g^ 

where /? = b^Dib. Since XDX ^ 0, we have 

/3|ci|^ - C > 0, (/3|cir - e)(/3 - ix") > (/3|ci|2 - ^xf (79) 
which results in |cip = With this, we can rewrite (l76l ) as 



cipbbt b 1 




^b^ 






bt X J 


141 




II 





.j-l II I 

which leads to X has rank one. This completes the proof. 
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Fig. 1. Eigenvalues of in iterations, positive definite case. 
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Fig. 2. Secrecy rate in iterations, positive definite case. 
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Fig. 3. ||Q*+^ - Q^ll, positive definite case. 
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Fig. 4. Secrecy capacity Vs. SNR, positive definite case. 
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Fig. 5. Eigenvalues of optimal input covariance Vs. SNR, positive definite case. 
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Fig. 6. Eigenvalues of in iterations, non positive definite case. 
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Fig. 7. Secrecy rate in iterations, non positive definite case. 
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Fig. 8. ||Q''+^ - Q*||, non positive definite case. 
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Fig. 9. Secrecy capacity Vs. SNR, non positive definite case. 
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Fig. 10. Eigenvalues of optimal input covariance Vs. SNR, non positive definite case. 
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